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Abstract 

Extending previous results of Oh-Zumbrun and Johnson-Zumbrun, we show that 
spectral stability implies linearized and nonlinear stability of spatially periodic traveling- 
wave solutions of viscous systems of conservation laws for systems of generic type, 
removing a restrictive assumption that wave speed be constant to first order along the 
manifold of nearby periodic solutions. 

1 Introduction 

Nonclassical viscous conservation laws arising in multiphase fluid and solid mechanics ex- 
hibit a rich variety of traveling wave phenomena, including homoclinic (pulse-type) and 
periodic solutions along with the standard heteroclinic (shock, or front-type) solutions 
[GZl \Z6\ IQZll IOZ2] . Here, we investigate stability of spatially periodic traveling waves: 
specifically, sufficient conditions for stability of the wave. 

In previous work [0Z44 I JZ3j . we showed that strong spectral stability in the sense of 
Schneider [SH \S2\ IS3| implies linearized and nonlinear L 1 n H — > L°° stability in all 
dimensions d > 1. However, as pointed out in [OZ1L ISel] . the conditions of Schneider are 
nongeneric in the conservation law setting, implying the restrictive condition that wave 
speed be constant to first order along the manifold of nearby periodic solutions. Indeed, 
it was shown in [OZ2] that failure of this condition implies a degradation in the decay 
rates of the Green function of the linearized equations about the periodic wave, suggesting 
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that nonlinear stability would be unlikely in the general (nonstationary wave speed) case 
in dimension d = 1. 

In this paper, we show that these difficulties are only apparent, and that, somewhat 
surprisingly, spectral stability implies nonlinear stability even if this additional condition 
on wave speeds is dropped. More precisely, we show that small L 1 n H s perturbations 
of a planar periodic solution u(x, t) = u(xi) (without loss of generality taken stationary) 
converge at Gaussian rate in L p , p > 2 to a modulation 

(1.1) u{ Xl - tP{x,t)) 

of the unperturbed wave, where x = (x\, x), x = (x2, ■ ■ ■ , Xd), and ip is a scalar function 
whose x- and t- gradients decay at Gaussian rate in all L p , p > 2, but which itself decays 
more slowly by a factor i 1 / 2 ; in particular, ip is merely bounded in L°° for dimension d = 1. 

In proving this result, we make crucial use of the tools developed in |OZ4j, IJZ3| . in 
particular, a key nonlinear cancellation argument of jJZ3j . The key new observation mak- 
ing possible the treatmen of the generic case is a careful study of the Bloch perturbation 
expansion about frequency £ = 0, motivated by relations to the Whitham averaged system 
observed in (Sell IJZTI IJZB] . 

It was shown in [Sel, OZ3| that the low- frequency dispersion relation near zero of the 
linearized operator about a periodic solution u agrees to first order with that of the lin- 
earization about a constant state of the Whitham averaged system 

d t M + Y,d Xj F J =0, 

(1-2) 

dtQiN) + v x (ns) = o, 

where M € M n denotes the average over one period, F 3 the average of an associated flux, 
Q = \V x y\ e M 1 the frequency, S = -^ t /\^x^\ G K 1 the speed s, and N = \7 X ^/\V X ^\ € 
M. d the normal v associated with nearby periodic waves, with an additional constraint 

(1.3) curl (UN) = curl = 0. 

As noted in [Sel|. IOZ3j , this implies both that the eigenvalues \j (£ ) bifurcating from A = 
at £ = are C 1 along rays through the origin, and that weak hyperbolicity (reality of 
characteristics of (jl.2|) — fjl.3|) ) is necessary for spectral or linearized stability. 

As noted in [i JZBj . there is a deeper analogy between the low- frequency linearized disper- 
sion relation and the Whitham averaged system at the structural level, suggesting a useful 
rescaling of the low-frequency perturbation problem. It is this intuition that motivates our 
derivation of sharp low-frequency estimates crucial to the analysis of nonlinear stability. 
With these estimates in place, the rest of the argument goes exactly as in | JZ3|. IOZ4| . 

1.1 Equations and assumptions 

Consider a parabolic system of conservation laws 

(1.4) ut + Y;fS(v) Xj =A x u, 

j 
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u G W(open) € R n , / J € M n , x E d > 1, t € M. + , and a periodic traveling wave solution 

(1.5) u = u(x ■ v — st), 

of period X, satisfying the traveling-wave ODE u" = - Ujf 3 (u))' — sv! with boundary 
conditions u(Q) = u{X) =: uq. Integrating, we obtain a first-order profile equation 

(1.6) u' = ^v j f :l (u) - su- q, 

j 

where (uq, q, s, v, X) = constant. Without loss of generality take u = ei, s = 0, so that 
u = u(xi) represents a stationary solution depending only on x\. 

Following [SeTllOZ3llOZ4] . we assume: 

(HI) p e C K+1 , K > [d/2] + 4. 

(H2) The map H : IxWxIx 5 d_1 x R n R n taking (X;a,s,u,q) ^ u(X;a,s,v,q)-a 
is full rank at (X;u(0),0,ei,q), where «(•; ■) is the solution operator of (ll.6p . 

Conditions (H1)-(H2) imply that the set of periodic solutions in the vicinity of u form 
a smooth (n + d+ l)-dimensional manifold {u a (x ■ is(a) — a — s(a)i)} , with a € M, a € R n+(i . 

1.1.1 Linearized equations 

Linearizing f)l .4[) about ti(-), we obtain 

(1.7) v t = Lv := A x v - ^(A> 



where coefficients A 3 := Df 3 (u) are now periodic functions of x\. Taking the Fourier 
transform in the transverse coordinate x, = (x2, ■ ■ ■ ,£d), we obtain 

(1. 8 ) v t = L^v = v XUXl - (A 1 v) Xl - i ^ZjV - Y £jV, 

where £ = (£2? • " >^d) is the transverse frequency vector. 

1.1.2 Bloch Fourier decomposition and stability conditions 

Following [Gl EH IS2| IS3J . we define the family of operators 

(1.9) L 5 = e~^ 1Xl L^ 1Xl 

operating on the class of I? periodic functions on [0, X}; the (L 2 ) spectrum of is equal 
to the union of the spectra of all with £1 real with associated eigenfunctions 

(1.10) w(xi,£, A) --e^qixu^^X), 

where q, periodic, is an eigenfunction of Lg. By continuity of spectrum, and discreteness of 
the spectrum of the elliptic operators on the compact domain [0, X], we have that the 
spectra of may be described as the union of countably many continuous surfaces Aj(£). 
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Without loss of generality taking X = 1, recall now the Bloch-Fourier representation 

(1.11) u(x) = (+-)* r [ e^u^x^di 

\2tt/ J^J^-i 

of an I? function u, where u(^,xi) := ^2 k e 2mkxi u(^i + 27r/c, £) are periodic functions of 
period X = 1, u(£) denoting with slight abuse of notation the Fourier transform of u in the 
full variable x. By Parseval's identity, the Bloch-Fourier transform u{x) — > u(£,xi) is an 
isometry in I?: 

(1-12) II u IIl 2 (x) = IHIl 2 (5;L 2 (:ei)), 

where L?(x\) is taken on [0,1] and L 2 (^) on [— tt,tt] x R^ -1 . Moreover, it diagonalizes the 
periodic-coefficient operator L, yielding the inverse Bloch-Fourier transform representation 

(1.13) e Lt u = (±-Y f / e^e^o^,^)^!^ 

V27T/ J-TrjRd-l 

relating behavior of the linearized system to that of the diagonal operators Lg. 

Loosely following [QZ4] . we assume along with (H1)-(H2) the strong spectral stability 
conditions: 

(Dl) a(L() C {ReA < 0} for f ^ 0. 

(D2) Rea(L ? ) < -6\£\ 2 : , 6 > 0, for £ G M d and |f | sufficiently small. 
(D3') A = is an eigenvalue of Lq of multiplicity exactly n + lQ 

As shown in [OZ3], (H1)-(H2) and (D1)-(D3') imply that there exist n + 1 smooth 
eigenvalues 

(1.14) A i (0 = -*a i (0 + o(iei) 

of bifurcating from A = at £ = 0, where — iaj are homogeneous degree one functions; 
see Lemma |2 . 1 1 below. 

As in [OZ4j, we make the further nondegeneracy hypothesis: 

(H3) The functions «-,(£) in (|1.14p are distinct. 
The functions a,- may be seen to be the characteristics associated with the Whitham av- 
eraged system (|1.2p - (|1.3jl linearized about the values of M, S, N, O, associated with the 
background wave u; see [OZ3| IOZ4] . Thus, (Dl) implies weak hyperbolicity of (|1.2p - (|1.3|) 
(reality of aj), while (HI) corresponds to strict hyperbolicity. 

Remark 1.1. Condition (D3') is a weakened version of the condition (D3) of [OZ4t |JZ3] 
that A = be a semisimple eigenvalue of Lq of minimal multiplicity n + 1, which implies 
[OZl[ IOZ2[ ISelj the special property that wave speed be stationary at u along the manifold 
of nearby periodic solutions. The stronger conditions (D1)-(D3) are exactly the spectral 
assumptions of |SH |S2"1 [S3] introduced by Schneider in the reaction-diffusion case. Condi- 
tions (D1)-(D3) (resp. (D1)-(D3')) correspond to "dissipativity" of the large-time behavior 
of the linearized system |Sll |S2"1 [S3] . 

1 The zero eigenspace of Lq is at least (n + l)-dimensional by linearized existence theory and (H2), and 
hence n + 1 is the minimal multiplicity; see [Sel| . 
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I. 2 Main result 

With these preliminaries, we can now state our main result. 

Theorem 1.1. Assuming (Hl)-(HS) and (D1)-(D3 7 ), for some C > andV> £ W K >°°(x,t), 

\u - u{- - rp)\u>(t) < C(l + tyU 1 - 1 ^ - u\ L ^ HK \t=u, 

(1.15) \u-u{- -tp)\ H K(t) < C(l + t)-i\u- u\ L i nH K\ t =o, 

\{&t,il>x)\wx+i* < C(l + tr% {1 ~ 1/p) \u - u\ L i nH K\ t=0 
for all t > 0, p > 2, d > 1, and 

(1.16) lu-Slip^), < C(l + t)"^ (1_ p )+ ^|n-n| L i nH K| t=0 

/or all t > and p = oo or p > 2 and d > 3, /or solutions u of (|1 .4[) [tt — u|^i n ^A- 1 t=o 
sufficiently small. In particular, u is nonlinearly bounded L 1 PI H — > L°° stable for d > I, 
asymptotically L 1 D — > L°° stable for d > 2, and asymptotically L 1 n fl" — > i? stable 
for d > 3. 

Remark 1.2. In Theorem \l.ll derivatives in x £ M d /or d > 2 re/er to total derivatives. 
Moreover, unless specified by an appropriate index, throughout this paper derivatives in 
spatial variable x will always refer to the total derivative of the function. 

In dimension one, Theorem 11.11 asserts only bounded L 1 n JET — > L°° stability, a very 
weak notion of stability. The bounds (|1.15p - (|1.16p agree for dimension d = 1 with those 
obtained in [JZ3] in the stationary wave speed case that (D3) holds in place of (D3'), but 
for higher dimensions are weaker by roughly factor t 1 / 2 . 

Remark 1.3. In dimension d = 1, it is straightforward to show that the results of Theorem 

II. 11 extend to all 1 < p < oo using the pointwise techniques of |OZ2j ; see Remark 13.61 



1.3 Discussion and open problems 

The proof of Theorem 11.11 largely completes the line of investigation carried out in [OZ2, 
ISeH IOZ3"l IOZ41 |JZ3| . showing that spectral stability implies linear and nonlinear stability of 
planar spatially periodic traveling waves. The corresponding spectral stability problem has 
been studied analyticially in |QZ1| ISeH IOZ3j . yielding various necessary conditions, and by 
numerical Evans function investigation in |OZl| . An interesting direction for further study 
would be more systematic numerical investigation along the lines of [BLZ, HLyZl, HLyZ2 



BHZl IBLZ] in the viscous shock wave case. A second interesting open problem would be 
to extend the results for planar waves to the case of solutions with multiple periods, as 
considered in the reaction-diffusion setting in |S1} IS2| [S3] . 

The key to the nonlinear analysis in critical dimensions d = 1, 2, as in |JZ3[ ISil IS2[ IS3|. 
is to subtract out a slower-decaying part of the solution described by an appropriate mod- 
ulation equation and show that the residual decays sufficiently rapidly to close a nonlinear 
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iteration. Note that the modulated approximation u{x\ — ip(x,t)) of (jl.ip is not the full 
Ansatz u a (^(x, t)), ^(x,t) := x\ — ip(x,t), associated with the Whitham averaged system 
(ll.2p - (ll.3p . where u a is the manifold of periodic solutions near u introduced below (H2), 
but only the translational part not involving perturbations a in the profile. (See [OZ3] for 
the derivation of Ansatz and (|1.2p -( fL3]l .) That is, we don't need to separate out all vari- 
ations along the manifold of periodic solutions, but only the special variations connected 
with translation invariance. 

This can be understood heuristically by the observation that (|1.2|) indicates that vari- 
ables a, are roughly comparable, which would suggest, by the diffusive behavior 
Vl/ >> V x ^/, that a is neglible with respect to Indeed, this heuristic argument translates 
rigorously to our ultimate computation of linearized behavior leading to the final result; 
see Section [2] and Remark I2.2I In this respect, the connection to the Whitham system is 
somewhat clearer in the generic case considered here than in the quasi-Hamiltonian case 
treated previously in [OZ2l IQZil IJZ3 ] Fl 

It would be interesting to better understand the connection between the Whitham av- 
eraged system (or suitable higher-order correction) and behavior at the nonlinear level, as 
explored at the linear level in |OZ3t lQZ4l IJZll |JZB] . As discussed further in [OZ3], another 
interesting problem would be to try to rigorously justify the WKB expansion for the related 
vanishing viscosity problem, in the spirit of [GMWZ1, GMWZ2J. 



2 Spectral preparation 

We begin by a careful study of the Bloch perturbation expansion at f = 0. 

Lemma 2.1. Assuming (H1)-(H3), (Dl)-(DS'), the eigenvalues Aj(f/|f|,f) of are an- 
alytic functions o/f/|f| and |f|. Suppose further that is a nonsemisimple eigenvalue of 
Lq, i.e., (D3 ! ) holds, but not (D3). Then, the Jordan structure of the zero eigenspace of 
Lq consists of an n- dimensional kernel and a single Jordan chain of height 2, where the 
left kernel of Lq is the n-dimensional subspace of constant functions, and u' spans the right 
eigendirection lying at the base of the Jordan chain. Moreover, for |f | sufficiently small, 
there exist right and left eigenf unctions <Zj(f/|f |, f , •) and <Zj(f/|f |, f , •) of associated with 
Xj of form qj = ^2 k (3j t kVk and qj = ^2 k /3j t kVk where {vj} and {vj} are dual bases of the to- 
tal eigenspace of associated with sufficiently small eigenvalues, analytic in cu = and 
|f |, with Vj(u;Q) constant for j ^ n and v n (u;;0) = u'(-); fy^, . . . , /3j >n -i, lll -1 ^, Pj, n +l 
and P jt i, ... ,/3j, n -i, If \Pj,n, Pj,n+l are analytic in f/|f|, |f|; and (qj,q k ) =5%. 

Proof. Recall that as an elliptic second-order operator on bounded domain has spectrum 
consisting of isolated eigenvalues of finite multiplicity. Expanding 

(2-1) Lt = L + \Z\Ll m + \Z\ 2 Ll m 



2 In the degenerate case that the stronger condition (D3) holds, i.e., wave speed is stationary at u, the 
situation is somewhat more complicated, and these relations break down; see [JZ3| for further discussion. 
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for each fixed angle £ := consider the continuous family of spectral perturbation 

problems in |£| indexed by angle to = about the eigenvalue A = of Lq. 

Because is an isolated eigenvalue of Lq, the associated total right and left eigenpro- 
jections Po and Po perturb analytically in both to and giving projection and P^ [K]. 
These yield in standard fashion (for example, by projecting appropriately chosen fixed sub- 
spaces) locally analytic right and left bases {vj} and {vj} of the associated total eigenspaces 
given by the range of P^, P^. 

Defining V = {v\, . . . ,v n+ i) and V = {v\, . . . ,v n+ \)* , * denoting adjoint, we may con- 
vert the infinite-dimensional perturbation problem (12. ip into an (n + 1) x (n + 1) matrix 
perturbation problem 

(2.2) M 5 = M + |£|Mi + |£| 2 M 2 , 

where M^(u, |£|) := (v£ , L^Vg) and (•, •) refers to the L 2 (xi) inner product on [0,X]. That 

is, the eigenvalues Aj(£) lying near of L^ are the eigenvalues of Mg, and the associated 
right and left eigenfunctions of L^ are 

(2.3) fj = Vwj and fj = wjV* , 

where Wj and Wj are the associated right and left eigenvectors of Mg. 

Case (i). If A = is a semisimple eigenvalue of Lq, then Mq = 0, and (|2.2p reduces 
to the simpler perturbation problem M ? := |^| _1 M ? = Mi + \t\M 2 studied in [OZ4[ IJZ3] . 
which = |£|Aj'(£)j Aj(^) denoting the eigenvalues of M^. Since Xj are continous, Aj 

are differentiable at |£| = in the parameter |^| as asserted in the introduction. Moreover, 
by (H3), the eigenvalues Aj (0) of Mi = Mq are distinct, and so they perturb analytically in 
to, |^|, as do the associated right and left eigenvectors. 

Case (ii). Hereafter, assume that A = is a nonsemisimple eigenvalue of Lq, so that 
Mq is nilpotent but nonzero, possessing a nontrivial associated Jordan chain. Moreover, as 
the n-dimensional subspace of constant functions by direct computation lie in the kernel of 
Lq = {d^ 1 +^1^), where Ai(x\) := df 1 (u(xi)), we have that the (n + l)-dimensional zero 
eigenspace of Lq is consists precisely of an n-dimensional kernel and a single Jordan chain 
of height two. Moreover, by translation-invariance (differentiate in x\ the profile equation 
(jl.6p ). we have Lqu' = 0, so that v! lies in the right kernel of Lq. 

Now, recall the assumption (H2) that H : IxWxIx S^ 1 x M n ->• W l taking 
(X; a, s, v, q) h-» u(X; a, s, v, q)—a is full rank at (X; u(0), 0, e±, q), where u(-; •) is the solution 
operator of profile ODE (]1.6p . The fact that kerLo is n-dimensional implies that the 
restriction H taking (a,q) i— > u(X;a, s,v,q) — a for fixed {X, v, s) is also full rank, i.e., 
H is full rank with respect to the specific parameters (X, s,v). Applying the Implicit 
Function Theorem and counting dimensions, we find that the set of periodic solutions, i.e., 
the inverse image of zero under map H local to u is a smooth (n + d + l)-dimensional 
manifold {u a (x ■ v(a) — a — s(a)t)}, with a G M, a E M n+rf . Moreover, d + 1 dimensions may 
be parametrized by (X, s, v), or without loss of generality (a\, . . . , a^+i) = (X, s, u). 
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Fixing (X, ^) and {ad+2-, ■ ■ ■ > a n+d+i)) and varying s, we find by differentiation of (jl.6p 
that := — 9 s u satisfies^! the generalized eigenfunction equation 

Lof* = u' . 

Thus, v! spans the eigendirection lying at the base of the Jordan chain, with the generalized 
zero-eigenfunction of Lq corresponding to variations in speed along the manifold of periodic 
solutions about u. Without loss of generality, therefore, we may take v%, . . . v n -\ and £„,+i 
to be constant at |£| = 0, i.e., depending only on u = £/|£| and not x\, and v n = v! at 
|£| = independent of u. 

Recalling from |JZ3| the fact that 

(c,LV) = (c, (^1(20^ - Ax) - J2"3 A j))u>') = {c^id 2 xi u -Y, u j d *J j (fi)) = 

for any constant functions c, where again (•, •} denotes L 2 {x\) inner product on the interval 
xi € [0, X], and Aj := we find under this normalization that (|2.2p has the special 

structure 

(2.4) M 
Now, rescaling (|2.2p as 

(2.5) M ? := ler^^M^^)- 1 , 
where 

(2.6) S :-- 
we obtain 

(2.7) M^Mo + lClM + O^ 2 ), 

where Mj = Mj(cv) like the original Mj are analytic matrix- valued functions of ui, and the 
eigenvalues rrij(£) = rrij(uj; |£|) of Mg are |£| -1 Aj(£). 

As the eigenvalues rrij of Mg are continuous, the eigenvalues Aj(£) = \£\mj are differen- 
tiable at |£| = as asserted in the introduction. Moreover, by (H3), the eigenvalues Aj(0) 
of Mo are distinct, and so they perturb analytically in oj, |£|, as do the associated right and 
left eigenvectors Zj and Zj. Undoing the rescaling (|2.5p . and recalling (|2.3p . we obtain the 
result. □ 

Remark 2.2. Note that the nth coordinate of vectors w € C n+1 in the perturbation 
problem (j2.2[) corresponds as the coefficient of u' to variations \& in displacement. Thus, 
rescaling (|2.5p amounts to substituting for \£ the variable ~ VE^ of the Whitham 

averaged system (11. 2j) . 




3 Note the function /» is X-periodic, and hence in the domain of Lo since we have fixed the period X. 
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3 Linearized stability estimates 

By standard spectral perturbation theory [K] . the total eigenprojection -P(£) onto the 
eigenspace of Lg associated with the eigenvalues Aj(£), j = 1, . . . ,n + 1 described in the 
previous section is well-defined and analytic in £ for £ sufficiently small, since these (by dis- 
creteness of the spectra of L^) are separated at £ = from the rest of the spectrum of Lq. 
Introducing a smooth cutoff function that is identically one for |£| < e and identically 
zero for |£| > 2e, e > sufficiently small, we split the solution operator S(t) := e Lt into 
low- and high-frequency parts 

(3.1) S J (tH := (^Y ^ [ e^mPiOe^MZ^iWidi 
and 

(3.2) S n (t)u := Q-) d £ J - ^P^e^uo^x^dl 



3.1 High-frequency bounds 

By standard sectorial bounds |Hel IPaj and spectral separation of Aj(£) from the remaining 
spectra of Lc, we have trivially the exponential decay bounds 

\\e L ^(I-cf>P(O)f\\ L m ,x]) < Ce- 0t \\f\\ L 2 {[o , x]) , 
(3-3) \\e L ^(I - <f>P(0)d l x Jh*(lo,x]) < Ct~ L 2e~ et \\f\\ LH[0)X]) , 

|lCe^(/-^P(C))/|| L2{M < Ct- L 2e- et \\f\\ L ,^ x]) , 

for 9, C > 0, and < m < K (K as in (HI)). Together with ()1.12p . these give immediately 
the following estimates. 

Proposition 3.1 QOZ4]). Under assumptions (H1)-(H3) and (D1)-(D2), for some 8, 
C > 0, and all t > 0, 2 < p < oo, < I < K + 1, < m < K , 

\\d l x s n {t)f\\ LHx): \\s n {t)d l x f\\L^) < ct-h- 9t \\f\\ LHx) , 

\\d™S n {t)f\\ LP{x) , W^^dTfh^ < Ct^-^e- et \\f\\ LHx) , 
where, again, derivatives in x <EM. d refers to total derivatives. 

Proof. The first inequalities follow immediately by f 1 1 . 1 2 1) . The second follows for p = oo, 
m = by Sobolev embedding from 

\\S n (t)f\\ Laaix . tLHxi)) < Ct-^e- 9t \\f\\ LH[0tX]) 

and 

\\d xl s II (t)f\\ L ^ {x . L 2 ixi)) < crV-i e - e *||/|| L2([0iX]) , 
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which follow by an application of (|1.12p in the x\ variable and the Hausdorff- Young in- 
equality ||/||i,ao(£) < ||/||x,i(D i n the variable x. The result for derivatives in x\ and general 
2 < p < oo then follows by LP interpolation. Finally, the result for derivatives in x follows 
from the inverse Fourier transform, equation (|3.2p . and the large |£| bound 

\e Lt f\m Xl ) < e~ mH \f\L2( Xl ), Ifl sufficiently large, 

which easily follows from Parseval and the fact that L^ is a relatively compact perturbation 
of d x — j£| 2 . Thus, by the above estimate we have 

\\e Lt dxf\\ LHx) < C\\e L ^\i\f\\ L 2 {xu0 

<Csup (e- e &\\) \\f\\ LHxi , 
<Ct^\\f\\ L2{x) . 

A similar argument applies for 1 < m < K. □ 

3.2 Low-frequency bounds 

Denote by 

(3.5) G I (x,t;y):=S I (t)5 y (x) 
the Green kernel associated with S 1 , and 

(3.6) [G|(xi,*;yi)] := mno^M] 

the corresponding kernel appearing within the B loch-Fourier representation of G , where 
the brackets on [Gg] and [S y ] denote the periodic extensions of these functions onto the 
whole line. Then, we have the following descriptions of G 1 , [Gg], deriving from the spectral 
expansion (|1.14p of L^ near £ = 0. 

Proposition 3.2 ([UZ4]). Under assumptions (Hl)-(HS) and (Dl)-(DS'), 

n+l 

[G|(xi, t; yi)] =<f>(0 £ e^&qjte, xjq^, yi)*, 

3=1 

(3.7) G\x,t ]y ) = (±-) d j S< x -^[Gl{ Xl ,t- yi )n 

where * denotes matrix adjoint, or complex conjugate transpose, qj(£, •) and qj(£, •) are right 
and left eigenf unctions of L^ associated with eigenvalues Aj(£) defined in (|1.14[) . normalized 
so that (qj ,qj) = 1 . 
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Proof. Relation (|3.T|) (i) is immediate from the spectral decomposition of elliptic operators 
on finite domains, and the fact that A,- are distinct for |£| > sufficiently small, by (H3). 
Substituting (13. 51) into (13.11) and computing 

(3.8) Syfox!) =J2 e2nikXl $y(£ + 2vrfce i) = ^2e 2nikxi e-^ y - 2nikyi = e-* v [5 yi (a*)], 

k k 

where the second and third equalities follow from the fact that the Fourier transform either 
continuous or discrete of the delta-function is unity, we obtain 

G\x,t-y) = (t^Y £ J /^(Oe^^xiR 

yielding (|3.7p (ii) by (|3.6p (i) and the fact that 4> is supported on [— tt, it]. □ 

Proposition 3.3. Under assumptions (H1)-(H3) and (D1)-(D3'J, the low-frequency Green 
function G I (x, t; y) of (|3.5p decomposes as G 1 = E + & , 

(3.9) E = u'(x)e{x,t;y), 
where, for some C > 0, all t > 0, 

supll^G,*,;^!^) <C(l + i)-* (1 -£> 

(3.10) V 

niT>\\%G J (;t,;v)\\i?(x), sup W t &{; t, ; y)\\ LP{x) < C{1 + t)-5< x -?H 
y y 

forp>2, 1 < r < 2, 

(3.11) S up||^^e(-,i,;y)|| iP(a;) < C7(l + t)-#Ci-|)-^-| 
/or p > 2, < j, k, I, j + 1 < K, 1 < r < 2, and 



(3.12) sup||ajd]e(-,i,;y)|| LP(:c) < C(l + i) 

y 



2^ r)' 2 



/or < j, k,l, j + I < K , provided that p > 2 and j + I > 1 or d > 3, or p = oo and d > 1. 
Moreover, e(x,t;y) = /or t < 1. 

Remark 3.4. In Proposition \3.3l and throughout the remainder of the paper, derivatives 
in y € M rf re/er to tota/ derivatives, just as with the variable x £ 
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Proof. In the degenerate case (D3) that is a semisimple eigenvalue of Lq, these estimates 
have been established in [0Z4t [JZ3]. Without loss of generality, therefore, we hereafter 
assume that is a nonsemisimple eigenvalue of Lq, with the consequences described in 
Lemma |2,1[ Recalling that 



G^x^y) = (±-Y [ e^yUiO^e^q^x^q^y,)*^ 

V27T/ J R d f— \ 

(3.13) 



1 \d 

2^ 



n+l 



define 



(3 . 14) e(x,t;y) = (-L) * / ^"^(O E ^MA^ViT^ 

R j,l 



so that 

G T (x,t;y) - u'(xi)e(x,t;y) 
1 \d 



(3.15) x2lT/ J^ d 



+ 



(i) V- 6<€ ' {X ~ V)m E eA ' ( ^i.^(^^^0-«X^i))«K^2/i)*^, 
where, by analyticity of w n , u n (^,a;i) — u'(xi) = 0(|£|), and so, by Lemma [2TT| 

(3.16) P j ,nP j ,i(v n (Z,x 1 ) -u'(xij)vi(S,K)* = 0(1) 
and 

(3.17) Pj,kPj,MZ>*iMZ>vi)* = for fc + n - 

Note further that v\ = constant unless Z = n, in which case {3ji = 0(|£|) by Lemma |2.1 
hence 



(3.18) ^ 1 (/3 J> /3 j7 ^ n (e,x 1 )-n'(x 1 ))^(e,yi)*j =0(|^|) 
and 

(3.19) ^(pj,kPjMZ>xiMt,ViY) = 0(\£\) for k + n. 

From representation (|3TT51) . bounds d3HSD d3HZD , and KAj (£) < -<9|£j 2 , we obtain by 
the triangle inequality 

(3.20) ||G J < C||e-*Kl a V(OIU*(0 <C(l + t)~l. 
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Derivative bounds follow similarly, since xi-derivatives falling on vjj- are harmless, whereas, 
by (|3.18p - (|3.19p . y\- or t-derivatives falling on vji or on e l ^'^ x ~ y > bring down a factor of |£| 
improving the decay rate by factor (1 + t) -1 / 2 . (Note that |£| is bounded because of the 
cutoff function (ft, so there is no singularity at t = 0.) 

To obtain bounds for p = 2, we note that (13.101 may be viewed itself as a Bloch- 
Fourier decomposition with respect to variable z := x — y, with y appearing as a parameter. 
Recalling (j!.12j) . we may thus estimate 
(3.21) 

sup||G 7 (-,t;y) - ue(-,t;y)\\ L 2 ix} < 



C ^ sup ||(/>(Oe Aj(€) *Wfc(^zi)^*(-,yi)?)i(-,yi)*||L2(5;L2( 2ie [ 0iX ] )) 
j,k^n,l y 

+ C > sup||0(f)e >™ [ n )vi{-,yi) Wl^^l^z^x])) 

< C ^ sup ||<K£) e ~ e|C|2t HL2 (e) sup ||ufc(-,zi)||L2 (0iX) ||^(-,yi)*|| LO o (0iX) 

+ C2^sup||0(€)e l?l ||l2 (5) sup||( — )\\L^{o,x)m{-,yi) ||l-(o,x) 

<C(1 + t)"f, 

where we have used in a crucial way the boundedness of v\ in and also the boundedness 
of 

{ j|j J ~^[V„(w;r) 

in L 2 , where < r < |£|. Derivative bounds follow similarly as above, noting that y- or t- 
derivatives bring down a factor |£|, while x-derivatives are harmless, to obtain an additional 
factor of (1 + i) -1 / 2 decay. Finally, bounds for 2 < p < oo follow by L p -interpolation. 
Defining 

(3.22) e{x,t;y) := X (t)e{x,t;y), 

where x 1S a smooth cutoff function such that x(t) = 1 for t > 2 and = for t < 1, 
and setting G := G — u'(xi)e(x,t;y), we readily obtain the estimates (|3.25p by combining 
the above estimates on G 1 — ue with bound (13.41) on C?^. 



4 This is clear for £ = 0, since Vj are linear combinations of genuine and generalized eigenfunctions, which 
are solutions of the homogeneous or inhomogeneous eigenvalue ODE. More generally, note that resolvent 
of Z/£ — 7 gains one derivative, hence the total eigenprojection, as a contour integral of the resolvent, does 
too- now, use the one-dimensional Sobolev inequality for periodic boundary conditions to bound the L°° 
difference from the mean by the (bounded) H 1 norm, then bound the mean by the L 1 norm, which is 
controlled by the L norm. 
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Finally, recalling, by Lemma 12.11 that vi = constant for I ^ n while f3j n = 0(|£|), we 
have 



Bounds (|3.1ip thus follow from (|3.14p by the argument used to prove (|3.10p . together with 
the observation that x- or i-derivatives bring down factors of |£|. 

Bounds (|3.12p follow similarly for p = oo if e~ e ^ t /\£,\ is integrable in ~R d , and for p > 2 
if e -6 *'^' Vl£| 2 i s integrable, thus yielding the stated results for all d > 2. In the special 
case d = 1, p = oo, (|2.1|) becomes a simpler one-parameter perturbation in £, and the 
contributions become analytic multiples of £ — whose principal value integrals may 
be carried out explicitly to give a sum of traveling error functions that is bounded in L°°; 
see the proof of Proposition 1.5, [OZ2] in the one-dimensional case. We omit this calculation 
as largely outside our analysis. (However, note that we need this bound to conclude L°° 
bounded stability in the one-dimensional case.) 

□ 

Remark 3.5. Underlying our analysis, and that of [OZ2[lJZ3] . is the fundamental relation 



(3.23) 



G(x,t;y) 



2irJ 



[G € (xi,t;yi)]d£. 



3.3 Final linearized bounds 

Corollary 3.1. Under assumptions (Hl)-(HS), (Dl)-(DS'), the Green function G(x,t;y) 
°f (|l-7p decomposes as G = E + G, 

(3.24) E = u'(x)e{x,t;y), 

where, for some C > 0, all t > 0, l<q<2<p< oo, < j, k, I, j + I < K , 1 < r < 2, 



(3.25) 



G(x,t;y)f(y)dy 



do 

+ 0O 



d r v G(x,t;y)f(y)dy 



+oo 



dZG(x,t;y)f(y)dy 



Lp(x) 
Lp(x) 

LP(x) 



xt-KV,-i/ 2 )-i| /k?ni2) 
<C(l + t)-i^ 2 ~ l M-\ +r 



xt-IWi-^-r\f\ LqnL2 . 



(3.26) 



oo 
oo 



d> x d«e(x,t;y)f(y)dy <(l+t) 



-|(i/ ff -i/ p )-fi+*i+l| / | iflj 



d> x $%e(x,t;y)f(y)dy 
Moreover, e(x, t; y) = for t < 1. 



LP 



< l+t)-2 



f(i/ g -i/ P )-(i±M 



I/Im- 
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Proof. (Case q = 1). From (|3. 10[) and the triangle inequality we obtain 

/ G I (x,t;y)f(y)dy <[ sup ||G'(., t; y)\\ LP \f(y)\dy < C(l + t)~i 

and similarly for y- and t-derivative estimates, which, together with (|3.4p . yield (|3,25p . 
Bounds (|3,26|) follow similarly by the triangle inequality and (|3.11j) — (|3.12p . 

(Case q = 2). From f|3.16p fj3.17|) . and analyticity of Vj, Vj, we have boundedness from 
L 2 [0, X] — > L 2 [0, X] of the projection-type operators 

(3.27) / -»• p jt J jtl (y n {£, xi) - u'(xi)) (5,, /) 
and 

(3.28) /^/3, )fe ^ fe (^xi)^,/)for k ^ n, 

uniformly with respect to £, from which we obtain by (|3.15p . (|3.22p . and f)1.12j) the bound 



(3-29) / G I (x,t;y)f(y)dy <C\\f\\ L2{x) , 

for all t > 0, yielding together with (|3.4|) the result (|3.25p for p = 2, r = 1. Similarly, by 
boundedness of Uj, i>j, in all L p [0, X], we have 

\e Xi(Ct MM^^iJ)\L^ Xl ) <Ce- e ^ 2t \f(i,-)\ L , {xi) , for k^n, 
C, 6 > 0, yielding by definitions (pTToT) . (I3T22]) the bound 

&(x,t;y)f{y)dy < (— ) / / C0(£)e~^ |2 'l/(£, •)Il=(«bi)«i ^ 

(3 ' 30) < c|^) e -^i 2 V (e3 l/lx-^o 

= C(l + t)-f||/|| x2([0|X])j 

hence giving the result for p = oo, r = 0. The result for r = and general 2 < p < oo then 
follows by LP interpolation between p = 2 and p = oo. Derivative bounds 1 < r < 2 follow 
by similar arguments, using (|3.18p - (|3.19p . Bounds (|3.26p follow similarly. 

(Case 1 < q < 2). By Riesz-Thorin interpolation between the cases q = 1 and g = 2, 
we obtain the bounds asserted in the general case 1 < q < 2, 2 < p < oo. □ 

Remark 3.6. The bounds on G, ej, e x may be recognized as the standard diffusive bounds 
satisfied for the heat equation [Z7j . For dimension d = 1, it may be shown using pointwise 
techniques as in |OZ2| that the bounds of Corollary 13.11 extend to all 1 < q < p < oo. 

We note a striking analogy between the Green function decomposition of Corollary 13.11 
and that of \MaZ3\ IZ4j in the viscous shock case; compare Proposition 3.3, [Z7| . 
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4 Nonlinear stability in dimension one 

With the bounds of Corollary (|3.ip , nonlinear stability follows by exactly the same argument 
as in |JZ3j . included here for completeness. We carry out the nonlinear stability analysis 
only in the most difficult, one-dimensional, case. The extension to the multi-dimensional 
case is straightforward |JZ3l IOZ4] . (Recall that the nonlinear iteration is easier to close in 
multi-dimensions, since the linearized behavior is faster decaying [OZ4t |JZ3l IS1L [S2l IS3|.) 
Hereafter, take x € R , dropping the indices on p and Xj and writing u t + f(u) x = u xx . 

4.1 Nonlinear perturbation equations 

Given a solution u(x,t) of (jl.4p . define the nonlinear perturbation variable 

(4.1) v = u — u = u(x + ij)(x, t)) — u(x), 
where 

(4.2) u(x, t) := u(x + tp(x, t)) 
and ip : R x R — > M. is to be chosen later. 

Lemma 4.1. For v, u as in (|4.ip . (|4,2|) . 

(4.3) ut + f{u) x - u xx = (d t - L) u'(xMx, t) + d x R + (d t + d 2 x )S, 
where 

1 2 I — I I I 

R : = wfo + V ip xx + (u x + v x )-p^- = 0(\v\(\ip t \ + \xfj xx \) + ( .I, ) IV'xf) 

1 + tPx \ 1- \lp x \ J 

and 

S:=-vip x = 0{\v\{\ip x \). 

Proof. To begin, notice from the definition of u in (|4.2p we have by a straightforward 
computation 

u t (x,t) = u x (x + if)(x,t),t)ip t (x,t) + u t (x + 
f(u(x,t)) x = df(u{x + ip(x,t),t))u x (x + ip,t) ■ (1 + i> x {x,t)) 

and 

^(ar,*) = (n x (x + if>(x,t),t) ■ (1 + if) x {x,t))) x 

= u xx (x + ip(x,t),t) ■ (l + ip x (x,t)) + (n x (x + ip(x,t),t) ■ ip x (x,t)) x . 

Using the fact that ut + df(u)u x — u xx = 0, it follows that 

^ 4 ^ u t + f(u) x - u xx = u x ip t + df(u)u x ip x - u xx ip x - {u x tp x ) x 

= U x Tp t ~ Utlpx ~ {U x ^ x ) x 
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where it is understood that derivatives of u appearing on the righthand side are evaluated 
at (x + ip(x,t),t). Moreover, by another direct calculation, using the fact that L(u'(x)) = 
by translation invar iance, we have 

(d t - L) u'(x)ip = u x tp t ~ u t ip x - (u x i/j x ) x . 

Subtracting, and using the facts that, by differentiation of (u + v)(x, t) = u(x + tj),t), 



(4.5) 
so that 

(4.6) 



u x + v x = u x (l + ip x ), 
u t + v t = u t + u x ip t , 



u x -u x -v x = -(u x + v x ) 



u t -u t -v t = -{u x + v x ) 



1 + ^' 



1 + 

we obtain 

ut + f(u) x - u xx = (d t - L)u'(x)t/j + v x ^ t ~ v t ip x - (v x ip x ) x + ((u x + v x ) ^ x ) , 

V 1 + lfj x J X 

yielding ([33]) by v x if> t - v t ip x = (vip t )x - (vip x ) t and (v x ifj x ) x = (vip x ) xx - (vip xx ) x . □ 
Corollary 4.2. The nonlinear residual v defined in (14. If) satisfies 

(4.7) v t -Lv = (d t - L) u'(x x )^ -Q x + R x + (d t + 8 2 X )S, 
where 

(4.8) Q := f(u(x + i,(x, t),t)) - f(u(x)) - df(u(x))v = 0(\v\ 2 ), 

(4.9) R := vip t + vip xx + (u x + 



1 + ll>x ' 

and 

(4.10) S :=-vip x = 0(\v\(\^ x \). 

Proof. Taylor expansion comparing fj4.3j) and ut + f(u) x — u xx = 0. □ 
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4.2 Cancellation estimate 

Our strategy in writing (|4.7|) is motivated by the following basic cancellation principle. 
Proposition 4.3 ( |HoZ] ) . For any f(y, s) G W n C 2 with f(y, 0) = 0, there holds 

rt 
10 



(4.11) J* J G(x,t - s;y)(d s - L y )f(y,s)dyds = f(x,t). 



Proof. Integrating the left hand side by parts, we obtain 
(4.12) 

G(x,0;y)f(y,t)dy - J G(x,t;y)f(y,0)dy + J j \d t - L y )*G(x,t - s;y)f(y, s)dy ds. 
Noting that, by duality, 

{d t - L y )*G{x, t- s;y) = S(x - y)5(t - s), 

5(-) here denoting the Dirac delta-distribution, we find that the third term on the righthand 
side vanishes in (|4.12|) . while, because G(x, 0;y) = 5(x — y), the first term is simply f(x, t). 
The second term vanishes by f(y,0) = 0. □ 

4.3 Nonlinear damping estimate 

Proposition 4.1. Let vq S H (K as in (HI)), and suppose that for < t < T, the H K 

norm of v and the H K (x,t) norms of tpt and ip x remain bounded by a sufficiently small 
constant. There are then constants 6\ t 2 > so that, for all < t <T, 

(4.13) \v(t)\ 2 HK < Ce- 0lt \v(O)\ 2 HK + cj\- e *- 8 ) (\v\ 2 L , + \(AM\ 2 H K [x>t) ) (s)ds. 

Proof. Subtracting from the equation (|4.4p for u the equation for u, we may write the 
nonlinear perturbation equation as 

(4.14) v t + {df{u)v) x - v xx = Q(v) x + u x ipt - u t ip x - (u x ip x ) x , 

where it is understood that derivatives of u appearing on the righthand side are evaluated 
at (x + ip(x, t),t). Using (|4.6p to replace u x and ut respectively by u x + v x — (u x + v x ) 1 ^ 

and tit + vt — (u x + v x ) j^ , and moving the resulting vttp x term to the lefthand side of 
(|4.14p . we obtain 



(1 + ip x )v t - v xx = -(df(u)v) x + Q(v) x + u x ip t 
{{ux + v x )ip x ) x + ((i 



+ \ {Ux+V x ) 



1 + V, 
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Taking the L 2 inner product in x of Ylf=o against (|4.15p . integrating by parts, and 
rearranging the resulting terms, we arrive at the inequality 

dt\v\ 2 HK {t) < -e\d« +1 v\ 2 L2 + c (\v\ 2 HK + mM\ 2 H K( X ,t)) > 

for some 6 > 0, C > 0, so long as \u\ h k remains bounded, and \v\ h k and \(tptj'4 , x)\H K (x,t) 
remain sufficiently small. Using the Sobolev interpolation \v\ 2 HK < Id^^v]^ + C\v\ 2 L2 for 
C > sufficiently large, we obtain d t \v\ 2 HK {t) < -6\v\ 2 HK + C (\v\\ 2 + \(t/j t , i>x)\ 2 H K( x ^ 
from which (14. 131) follows by Gronwall's inequality. □ 

4.4 Integral representation/ ^-evolution scheme 

By Proposition H31 we have, applying Duhamel's principle to (|4.7|) . 

poo 

v(x,t) = I G(x,t;y)v (y)dy 

t roo 



(4.16) ■ x 



+ 



ft f'OO 

/ / G(x,t-s;y)(-Qy + R x + St + Syy)(y,s)dyds + ip(t)u'(x). 

JO J-oo 



Defining ip implicitly as 
(4.17) 



oo 

ip(x,t) = -l e(x,t;y)u (y)dy 

oo 

t r+oo 



ft p+oo 

/ / e(x,t-s;y)(-Q y + R x + S t + Syy)(y,s)dyds, 

JO J-OD 



following [ZHl [Z4l [MaZ2l IMaZ3j . where e is defined as in (pT24"D . and substituting in ([4TT6]) 
the decomposition G = u'(x)e + G of Corollary 13. 1\ we obtain the integral representation 



oo 



v(x,t) = I G(x,t;y)v (y)dy 
t roo 



(4.18) 



ft roo 

+ G(x,t- s;y)(-Q y + R x + S t + S yy )(y,s)dyds, 

Jo J -oo 

and, differentiating (|4.17p with respect to t, and recalling that e(x,s;y) = for s < 1, 



%d*i/)(x,t) = - I d 1 t d^e(x,t;y)u (y)dy 

(4.19) 



oo 

ft f + OO 

/ / d{d^e(x,t- s;y)(-Q y + R x + S t + S yy )(y,s)dyds. 

Jo J-oo 



Equations (|4.18p . (|4.19j) together form a complete system in the variables (v, dfip, d x ip), 
0<j<l,0<k<K, from the solution of which we may afterward recover the shift ip via 
(|4.17p . From the original differential equation (|4,7p together with (|4.19p . we readily obtain 
short-time existence and continuity with respect to t of solutions (v^ti^Px) £ H by a 
standard contraction- mapping argument based on (|4.13j) . (|4.17j) . and and (|3.26p . 
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4.5 Nonlinear iteration 

Associated with the solution (u,^ t ,ip x ) of integral system (|4,18|) - (|4.19|) . define 
(4 20) ({*)'■= SU P \{ v ^t^x)\HK{s){\ + s) 1/4 . 

Lemma 4.2. For all t > for which ((t) is finite, some C > 0, and Eq := \uo\ L i nH K , 

(4.21) at)<c(E +at) 2 ). 

Proof. By (l49j) - (l4T0ll and definition (14201) . 

(4.22) \(Q,R,S)\ L i nL oo < ^, Vt, + \{v,v x ^ t ,tp x )\ 2 Laa < C((t) 2 (l + t)~^ , 

so long as \if) x \ < \ip x \ H K < remains small, and likewise (using the equation to bound t 
derivatives in terms of x-derivatives of up to two orders) 

(4.23) \(d t + dl)S\ L i nL oo < \(v^ x )\ 2 H2 + 1(^,^)1^,00 <CC(t) 2 (l + t)-i 

Applying Corollary 13.11 with q = 1, d = 1 to representations (|4.18p - ([4.19p . we obtain for 
any 2 < p < oo 

H-,t)\ LPix) <c(i+t)-^-yp) Eo 

+ c((t) 2 [ (i + t-sy^-^it-sy-^i + syUs 

Jo 

icciEo+atf^i+t)-^ 1 - 1 ^ 



(4.24) 



and 
(4.25) 

\{thM(;t)\ WKtP < c{\ + t)-7Eo + cat) 2 [ (i + 1 - s )-K 1 - 1 /p)-i/2 (1 + s) -| ds 

JO 

iCCiEo + atf^l + t)-^ 1 - 1 ^. 

Using (I4TT3D and (Oil) - K25\\ . we obtain t)\ H K^ < C(E + ((t) 2 )(l+t)-^ . Combining 
this with (|4.25p . p = 2, rearranging, and recalling definition (|4.20|) . we obtain (|4.2|) . □ 

Proof of Theorem By short-time iJ existence theory, ^(v^t^x 

)\\ H K is continuous 

so long as it remains small, hence r\ remains continuous so long as it remains small. By 
(|4.2p . therefore, it follows by continuous induction that n(t) < 2Ct]q for t > 0, if rjo < l/4:C, 
yielding by (|4T20|) the result (fLl"5|) for p = 2. Applying ()4T24"l) - (l435l) . we obtain (fTTT5l) for 
2 < p < for any < oo, with uniform constant C. Taking p* > 4 and estimating 

|Q| L2 , |i?| i2 , |S| L2 (i) < |(«,^,^)|i4<C^(l + t)-t 
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in place of the weaker (|4.22p . then applying Corollary 13.11 with q = 2, d = 1, we obtain 
finally (|1.15p for 2 < p < oo, by a computation similar ()4.24j) — ()4.25|) ; we omit the details of 
this final bootstrap argument. Estimate (I1.16P then follows using (13.260 with q = d = 1, by 

\m\LP < ce (i + t)a? + cc(o 2 / (i +t - sy^-yp)^ + s )-3ds 

(4.26) Jo 

<C(l + t)i(S + CW 2 ), 

together with the fact that u(x, t) — u(x) = v(x — ip,t) + (u(x) — u(x — ip), so that |n(-, t) — u\ 
is controlled by the sum of \v\ and \u(x) — u(x — ip)\ ~ This yields stability for 

\u — u\ L i nH K \t=o sufficiently small, as described in the final line of the theorem. □ 
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